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ABSTRACT 

Measurements of rates of period change of Classical Cepheids probe stellar physics and evolution. 
Additionally, better understanding of Cepheid structure and evolution provides greater insight into 
their use as standard candles and tools for measuring the Hubble constant. Our recent study of 
the period change of the nearest Cepheid, Polaris, suggested that it is undergoing enhanced mass 
loss when compared to canonical stellar evolution model predictions. In this work, we expand the 
analysis to rates of period change measured for about 200 Galactic Cepheids and compare them to 
population synthesis models of Cepheids including convective core overshooting and enhanced mass 
loss. Rates of period change predicted from stellar evolution models without mass loss do not agree 
with observed rates whereas including enhanced mass loss yields predicted rates in better agreement 
with observations. This is the first evidence that enhanced mass loss as suggested previously for 
Polaris and 6 Cephei must be a ubiquitous property of Classical Cepheids. 
Subject headings: stars: evolution — stars: variables: Cepheids — stars: mass-loss 



1. INTRODUCTION 

Observations of real-time stellar evolution are rare and 
difficult to obtain as most stars change on time scales 
much longer than a human life time. However, stel- 
lar pulsation proper ties are observed to change. Just 
as lEddingtonl ()1926D demonstrated that observations of 
stellar pulsation and the pulsation period constrain the 
interior structure of a star via the period-mean den- 
sity relation, changes in the pulsation period constrain 
the evolution of that star as changes in the mean den- 
sity. Currently, period change has been measured for 
about 200 Classical Cepheids ()Turner e t al.' '20 06D as 
well as a number of Type II Cep heids (Schmid t et al.l 
12004', '2005; 'Rabidoux ct al.' '201 0|) and RR Lyra e stars 
(|Le Borgne et al. 2007: Vandenbroere et al.l[2012l ). 

Rates of period change constrain which stage a Cepheid 
is evolving: a negative rate of period change suggests a 
Cepheid is evolving blueward on the Hertzsprung- Russell 
diagram on its second crossing of the Cepheid instabil- 
ity strip, while a positive rate suggests redward evolu- 
tion on the first or third crossing. The first crossing 
occurs when a star evolves through the Hertzsprung gap 
at the end of main-sequence evolution. A star cools as 
it crosses the Hertzsprung gap on a thermal time scale 
hence crosses the Cepheid instability strip more quickly 
than in later stellar evolutionary phases which are on a 
nuclear timescale. The rate of period change of a Cepheid 
in this phase of evolution is positive and greater than 
rates of period change during later crossings of the in- 
stability strip. The third crossing occurs near the end of 
helium-core burning when the star evolves to the asymp- 
t otic giant branch or t o the red supergiant branch. 

iTurner et al.l (|2006( ) found that measured pulsation 
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periods and rates of period change appear consistent 
with predic tions from stellar evolution models. However, 
iNeilson et~ al. (2012) computed a grid of stellar evolution 
models to fit the observed effective temperature and lu- 
minosity of Polaris but could not simultaneously fit its 
observed rate of period change. Predicted rates of period 
change consistent with the observed luminosity and ef- 
fective temperature are too small by more than 1 s yr^^. 
Stellar mass loss at the rate of M ~ 10~^ Mq yr~^ 
was suggested to account for the difi'erence between the- 
oretical and observed period changes. This hypothesis 
is contentious, as the necessary mass loss rate is rather 
high, but if correct is an important key to understanding 
Cepheid structure and evolution. 

Cepheid mass loss is an important ingredient for un- 
derstanding late-stage stella r evolution, the tra nsitions 
from red-to-blue supergiants (|Mackev et al.ll20l"2l ) and to 
solve the l ong-s tandi ng problem of the Cep heid mass dis- 
creoa ncv (iCoxl 11980^ 'Brocato et al.' '2004'; 'Caputo et alj 
2005: lBono eta l. 2006; Keller 2008; Ncilson ct al. 2011). 
The Cepheid mass discrepancy is the difference between 
Cepheid mass predictions using stellar evolution mod- 
els and stellar pulsation models, where the latter pre- 
dicts masses approximately 10-20% smaller. iBono et al.l 
(2006) suggested a number of potential solutions to re- 
solve the mass discrepancy with the two most promising 
ones being Cepheid mass loss and convective core over- 
shooting in main-sequence Cepheid progenitors. Mass 
loss removes mass from the stellar envelope whereas over- 
shooting leads to a more luminous Cepheid for a given 
i nitial st ellar mass. 

i Keller l (12008) argued that the mass discrepancy in- 
creases for decreasing metallicity, which is inconsistent 
with radiative-driven mass loss theory, hence advocat- 
ing convective c ore overshooting a s the s olutio n. On 
the other hand, INeilson fc Lesteii (|2008l I2009D devel- 
oped a pulsation-driven mass loss model and found that 
pulsation-driven mass loss does not decrease with de- 
creasing metallicity. Neilson et al.l (|2011[ ) added the 
pulsation-driven mass-loss theory to the stellar evolution 
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models and found that a combination of overshooting 
and Cepheid mass loss is required to account for the mass 
discrepancy. 

Infrared and radio observations also provide evidence 
for Cepheid mass loss. Interferometric observations 
of a number o f Cepheids suKgest if-band flux excess 
(jKervella et al.ll2006l: iMerand et al.ll2Q06ll20Q7[ ). indicat- 
ing the presence of circumstellar material from a stel- 
lar wind, but no flux excess about the non-pulsating 
yellow supergiant a Persei. Infrared excess was also 
found in iggztzer observations of Galac tic Cepheids 
(|Marengo et 'all 1213 ioi : fBarmbv et al] fMll ). in particu- 
lar, observations resolv e a bow shock about t he Ce pheid 
prototype & Cephei. iNeilson et all (|2009l I2O1O0 also 
found evidence for infrared excess due to mass loss 
in Large Mage ll anic Cloud Cepheids. Furthermore, 
iMatthews et al.l ()2012l ) detected significant mass loss in 
5 Cephei based on 21 cm radio observations, includ- 
ing a measurement of the wind velocity. All of these 
analyses suggest mass-loss rates of the order 10~^ - 
10^^ Mq yr~^, but could be explained by other causes 
such as material left over from other stages of evolu- 
tion. These mass-loss rates are much larger than is ex- 
pecte d from a radiatively-dr iven wind or a Reimers rela- 
tion (jNeilson fc Lesteiil2008D . 

Measurements of the period change of Polaris provided 
a first direct test of Cepheid mass loss. Whereas ob- 
servational evidence for Cepheid mass loss is based on 
flux excess, the analysis of the rate of period change 
measures how mass loss changes the evolution of the 
Cepheid. Period change also varies differently as a func- 
tion of convective core overshoot and mass loss, providing 
the first method t o constrain differen t solutions to the 
mass discrepancy (jNeilson et al.l [201^ . However, mea- 
suring mass loss from period change for Polaris is not 
reliable because Polaris is only one star and its large rate 
of period change may be explained by o ther phenomena, 
such as magnetic fields (|Stothersll2009D . If mass loss is 
so significant then there must be some evidence for it 
in the measured period changes of all Galactic Cepheids 
(lSzabados"1983': 'T urner et al.l l2006D. 

INeilson fc Lesteri (|2008| ) and INeilson et al.l (j2012l ) 
showed that the rate of period change is linearly pro- 
portional to the mass-loss rate, i.e. if mass loss increases 
then so too must the rate of period change. If Cepheid 
mass loss is ubiquitous then one would expect stellar evo- 
lution models to underestimate rates of period change 
relative to the measured rates as well as to predict more 
Cepheids wit h a negative rate of period change than a 
positive rate. I Turner et al.l (|2006D did not find evidence 
for the first hypothesis, however they did not conduct 
a quantitative comparison. The purpose of this arti- 
cle is to compare the me asured rates of period changes 
from I Turner et al.l (|2006[ ) to predicted rates computed 
from population synthesis models of Galactic Cepheids 
to test whether evolution models are consistent with the 
observed ratio of Cepheids with positive and negative pe- 
riod change. In Sect. 2, we describe the method for com- 
puting population synthesis models and the underlying 
model assumptions. In Sect. 3 we compare the measured 
rates of period change with stellar evolution models as- 
suming no significant Cepheid mass loss and repeat the 
comparison in Sect. 4 but for evolution models with sig- 
nificant Cepheid mass loss. We discuss the implications 



of this work and conclude in Sect. 5. 

2. METHODS 

\ ye compute detailed s tellar evolution models using 
the'Yoon & Langer (j2005l ) code, as in our previous work 
(iNeilson et al., 2012j). Here, we use that grid of models 
with masses, M = 3.2 - 6.8 M©, in steps of 0.1 Mq, with 
convective core overshooting para meter etc = - 0.4 in 
steps of 0.1 and assuming the Greves se et al.l (|1996D stan- 
dard solar metallicity. Convective core overshooting is 
discussed in Neilson_el al. (2011, 201^. We extend this 
grid up to 14 Mq . Rates of period change are computed 
from the stellar evolution tracks when the star's effective 
temperature and luminosity a re consistent with crossing 
the Cepheid instability strip (jBono "eraI|[2Q00 f) and for 
luminosities, lo^L/Lp, > 3, where rates a re computed 
using Equation 2 from INeilson et all (|2012l ). 
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where M is the change of stellar mass. 

We then compute population synthes i s mod els for 
Galactic Cepheids, assuming a iKroupal (|200lD initial 
mass function and a constant star formation rate. From 
our models, we compute the relative probability of see- 
ing a Cepheid with luminosity within a bin L + dL and 
L — dL and period change less than P + dP and greater 
than P — dP. By summing all the Cepheids with positive 
period change and those with negative period change, we 
compute the relative fraction of each that can be com- 
pared to observations. 

This comparison allows us to test the hypothesis that 
Cepheids undergo enhanced mass loss. Because mass loss 
increases the rate of period change then one would intu- 
itively expect population synthesis models would predict 
that the fraction of Cepheids with positive period change 
to be less than the observed fraction from lTurner et ahl 
(|2006f) . However, mass loss may also change evolutionary 
time scales, which would a lso ch ange the relative num- 
bers. In the iTurner et al.l (|2006( ) sample, there are 196 
Galactic Cepheids with measured rates of period change, 
of which 128 Cepheids are measured to have positive 
rates of period change and 68 negative rates of period 
change. Therefore, approximately 2/3 of all Cepheids 
have a positive rate of change. 

3. POPULATION SYNTHESIS MODELS WITHOUT 
CEPHEID MASS LOSS 

We compute rates of period change for our grids of stel- 
lar evolution models. The largest positive rates of period 
change are for stars on the first crossing during which the 
crossing time scale is proportional the the thermal time 
scale of the star. The second and third crossings take 
much longer, which lead to smaller absolute values for 
the rate of period change. We show the time scales of 
evolution for a 5 Mq model in Fig. |l] in which we plot 
the Hertzsprung-Russell diagram of the track and plot 
points at every time step of 10^ yr. The star evolves along 
the first crossing in about 1 Myr and the third crossing 
appears to be much longer than the second crossing, in 
apparent contradiction to our intuitive hypothesis that 
there would be fewer Cepheids with positive rates of pe- 
riod change, i.e. third-crossing Cepheids. 
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Figure 1. The Hertzsprung-Russell diagram for a 5 Mq stellar 
evolution track from our sample where the points represent time 
steps of 10^ yrs of stellar evolution. Dott ed lines re present the 
boundaries of the Cepheid instability strip IIBono et aI.ll200Q ). 

We explore this result further by predicting the proba- 
bility for observing a Cepheid with a positive or a nega- 
tive rate of period change as a function of luminosity from 
stellar evolution models without convective core over- 
shooting and plot the results in Fig. [2] There are three 
apparent tracks in the two plots representing the first and 
third crossings of the instability strips in the left plot and 
the second crossing in the right plot. Measured rates o f 
period change for Galactic Ce pheids (iTurner et al. 120061 ) 
are plotted for comparison. iTurner et alJ (|2d06 ) pre- 
sented the measurements as P = f{P), which we use 
in this work. However, we convert the pulsation periods 
to bol ometric luminosities using the relation from lTurneii 
((Moh for more direct comparisons with stellar evolution 

models as P/P = f{L/LQ). 

The observed rates of period change do not agree well 
with predicted rates for a given luminosity: predicted 
positive rates tend to be larger than the observed rates, 
especially for luminosities log L/Lq > 3.5 while observed 
negative rates have more scatter than predicted rates. 
This difference is resolved if we employ stellar evolution 
models that assume a moderate amount of convective 
core overshoot, ac — 0.2. Convective core overshoot acts 
to shift the contours to greater luminosities but main- 
taining the same behavior. We plot the rates of period 
change for the ac — 0.2 stellar evolution models in Fig.[3l 
The predicted positive rates of period change for mod- 
els with overshoot agree better with measured rates for 
logL/L© > 3.5, but disagree with period changes for 
Cepheids with smaller luminosities. Including convec- 
tive core overshoot in stellar evolution models increases 
a Cepheid's luminosity but the rates of period change do 
not appear to change. 

We calculate the fraction of Cepheids that have posi- 
tive period change and the fraction with negative period 
change from the synthesis models. For the stellar evo- 
lution models without overshooting, we find that about 
85% of Cepheids have positive rates of period change. 
We find a similar fraction of Cepheids from the stellar 
evolution models that include convective core overshoot. 
These two results differ from the observed fraction of 
Cepheids with positive rates of period change, about 67% 
and suggest that convective core overshoot does not af- 
fect the time scale for a Cepheid evolving redward rela- 



tive to its time scale for the second crossing. A different 
physical process is needed. 

While it is expected that Cepheid mass loss increases 
the relative fraction of Cepheids with positive rates of 
period change, this assumes that mass loss, like convec- 
tive core overshoot, does not affect the relative evolu- 
tionary time scales of the second and third crossings for 
a given mass. Previous tests suggest that mass loss does, 
in fact, change the relative time scales as well as the 
shape of t he blue loop on the Hertzsprung-Russell dia- 
gram (Ne ilson et al.l 120111 : iMatthcws et al. 2012). It is 
necessary to test whether Cepheid mass loss increases or 
decreases the predicted fractions of Cepheids with posi- 
tive and negative rates of period change by changing the 
relative evolutionary time scales of blueward and redward 
evolution. 

4. POPULATION SYNTHESIS MODELS WITH ENHANCED 
CEPHEID MASS LOSS 

Convective core overshoot does not resolve the dis- 
crepancy between the predicted and observed ratios of 
the number of Cepheids with positive rates of period 
change and negative rates of period change. We repeat 
the analysis from the previous section with a grid of stel- 
lar evolution models computed assuming no convective 
core overshoot while including enhanced mass loss at a 
rate M = 10"'' Mq yr""'^ when the stars evolve across 
the Cepheid instability strip. 

We show in Fig. |4] the rates of period change plotted 
as a function of stellar luminosity for the stellar evo- 
lution models with enhanced Cepheid mass loss. Pe- 
riod change predictions for the second crossing do not 
change significantly (Fig. HI right panel) because of the 
enhanced mass loss. However, predicted positive rates of 
period change differ, especially rates of period change of 
Cepheids evolving along the third crossing of the instabil- 
ity strip. Cepheid models evolving on the third crossing 
with luminosities log L/Lq = 3.2 - 3.5 have rates of pe- 
riod change up to one to two orders-of-magnitude greater 
than those predicted from models assuming no enhanced 
mass loss. The differences at larger luminosities may be 
resolved by including convective core overshooting in the 
stellar evolution calculations, as shown in Fig.O As mass 
loss appears to increase the rate of period change then 
massive Cepheids could have smaller mass-loss rates. 

From the models including enhanced mass loss, 71% 
of the Cepheids have positive rates of period change, a 
significantly smaller fraction than that predicted from 
stellar evolution models without mass loss enhancement, 
85%, and in better agreement with that measured from 
observations, 67%. This result suggests that theoretical 
stellar evolution models require enhanced mass loss on 
the Cepheid instability strip to predict rates of period 
change consistent with observations. 

5. DISCUSSION 

We computed grids of stellar evolution models for three 
cases: 1) with no convective core overshooting and no 
enhanced Cepheid mass loss, 2) convective core over- 
shooting, ttc = 0.2 and no enhanced mass loss and 3) 
no overshooting but with enhanced Cepheid mass loss, 
M = 10~^ Mq yr~^. For each grid, population synthe- 
sis models are used to predict the fraction of Cepheids 
with positive rates of period change and those with nega- 
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Figure 4. Same as Fig. [2] except using stellar evolution models with mass-loss rates M = 10 Mq yr ^ while evolving on the Cepheid 
instability strip. 



tive rates of period change. These predicted fractions are 
compared to the fraction of Cepheids with positive and 
negat ive rates measured from observations () Turner et aLl 
l2006l) . 



In the first two cases our models overestimate the frac- 
tion of Cepheids with positive rates of period change and 
convective core overshooting does not appear to change 
the fraction. However, in the third case, the fraction of 
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Table 1 

Fractions of Cepheids with positive and negative rates of period 
change 



Case 


Positive 


Negative 




85% 


15% 


€ 


85% 


15% 


3^ 


71% 


29% 


Observed 


67% 


33% 



''Models without overshoot or enhanced mass loss. 
''Models with overshoot but no enhanced mass loss. 
•^Models with enhanced mass loss but no overshoot. 

Cepheids with positive period change is smaller, improv- 
ing agreement with the measured fraction. We summa- 
rize the predicted ratios in Tab. [TJ This result is the 
first direct suggestion that Cepheids, in general, must 
undergo enhanced mass loss. Mass loss has been inferred 
from excess flux in the infrared and radio observations, 
but these could also be interpreted as a relic from previ- 
ous epochs of stellar evolution. The velocity of the cir- 
cumstellar medium of 6 Cephei was also measured from 
radio observations, suggestive of but not a direct mea- 
surement of a Cepheid wind. 

The strongest evidence for Cepheid mass loss would 
be the discovery of a P Cygni profile from spectral ob- 
servations. P Cygni line profiles may not be detectable 
in Cepheids as the spectral lines arc dominated by pulsa- 
tion and turbulence (Bcrsicr & Burki 1996). Therefore, 
rates of period change appear to provide the most direct 
evidence for enhanced mass loss in Cepheids. 

In this analysis we assumed a mass-loss rate M = 
10~^ yr~^, rou ghly consistent with t he predictions 
from iNeilson fc Les ter' (2008, 2009) and iNeilson et all 
(|2011t ). but without physical justification. This partic- 
ular average mass-loss rate leads to a more significant 
increase of the rate of period change in dimmer, i.e less 
massive M < 7 Mq, Cepheids but the change is negli- 
gible at high er mass. We c a ll this the average mass-loss 
rate because INeilson et al.l (|2011[ ) suggested that mass 
loss may change by orders of magnitude as a Cepheid 
evolves across the instability strip. There are fewer mas- 
sive Cepheids than smaller mass Cepheids because of the 
shape of the initial mass function, hence the smaller num- 
ber of massive Cepheids contribute little to the result- 
ing fraction of Cepheids with positive or negative period 
change. 

This mass loss is consistent with measu r ement s of 
the Cepheid mass discrepancy by iKelled (|2008l ). if 
we were also include some convective core overshoot- 
ing. For instance, measurements of the Large Mag- 
ellanic Cloud ecl ipsing binary OGLE-LMC-CEP0226 
(|Pietrzvhski et aLl|2010i) suggest that the Cepheid com- 
pone nt has convective core overshooting with a,. = 
0.2 (ICassisi fc Sakrisl I20TTI : INeilson fc La^^ I20T2I : 
iPrada Moroni et al.ll20"m 

While mass-loss rates of 10""^ Mq yr~^ are consis- 
tent with observations, it is unclear how a Cepheid loses 
mass at such high rates. Radiative- driven wind the- 
ory p redicts rates < 10~^° Mq yr~^ (jNeilson fc Lesteii 
I2OO8I ) yet Cepheids are too hot for dust-driven winds to 
form. INeilson fc Lesteii (|2008l I2OO90 suggested an ana- 
lytic model for pulsation-driven mass loss in Cepheids 



but the predicted mass-loss rates may still be too small 
to fit more recent observations, in particu lar 6 Cephei 
(jMarengo et al.ll2010HMatthews et al.ll2012D . A new de- 
scription for Cepheid mass loss is required to better agree 
with mass-loss rates inferred from observations as well as 
to better understand measured rates of period change. 

We demonstrate that mass loss is important for 
Cepheid evolution and thus might also affect evolution 
along the asymptotic giant branch. While our analysis 
is primarily applicable to intermediate mass stars, sig- 
nificant mass loss in massive stars, M > 8 Mq, might 
lead to lower mass supernova progenit ors (jEldridge et al.l 
l2008tlMaund et al.l201lllGeorgvl20il) . and could change 
the blue-to-red supergi ant ratio that is not understood 
from theoretical models (jMaeder fc Mevn"etll2012l : lLangeii 
l2012f l. Along with impacting stellar evolution, Cepheid 
mass los s may also affect the Cepheid period-luminosity 
relation. INeilson et a l. (2010) suggested that infrared ex- 
cess caused by mass loss can change the relation by 1-2%. 
Understanding the physical mechanism for mass loss in 
Cepheids is important for constraining the distance scale 
and for understanding stellar evolution of intermediate 
and massive stars. 
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